Effective mass in cavity QED 
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We consider propagation of a two-level atom coupled to one electro-magnetic mode of a high-Q 
cavity. The atomic center-of-mass motion is treated quantum mechanically and we use a standing 
wave shape for the mode. The periodicity of the Hamiltonian leads to a spectrum consisting of 
bands and gaps, which is studied from a Floquet point of view. Based on the band theory we 
introduce a set of effective mass parameters that approximately describe the effect of the cavity on 
the atomic motion, with the emphasis on one associated with the group velocity and on another one 
that coincides with the conventional effective mass. Propagation of initially Gaussian wave packets 
is also studied using numerical simulations and the mass parameters extracted thereof are compared 
with those predicted by the Floquet theory. Scattering and transmission of the wave packet against 
the cavity are further analyzed, and the constraints for the effective mass approach to be valid are 
discussed in detail. 



PACS numbers: 45.50.Dv, 42.50.Pq, 42.50.Mk 

I. INTRODUCTION 

Cavity quantum electro dynamics (QED) has expe- 
rienced a tremendous progress during the last decades. 
In experiments where atoms interact with a cavity field, 
the lifetimes of both the cavity and atomic states can be 
made rather long, up to tens of milliseconds. This makes 
it possible to carry out several operations on the com- 
bined system before decoherence plays an influential role. 
It is also possible to single out a unique atomic transition 
to interact with only one cavity mode, implying that only 
two atomic states | j) and | J.) together with one electro- 
magnetic mode \n) need to be taken into account, while 
other atomic states and modes are neglected. In such 
situations, the Jaynes-Cummings (JC) model 0, 0] has 
proven to be remarkably well suited for describing the 
coherent interaction. Cavity QED has thus become one 
of the candidates for implementing quantum information 
processing, see for example 0, U IE an d h has also 
turned out to be a very useful tool for studying purely 
quantum mechanical phenomena, such as sub-poissonian 
Fock states || and Schrodinger cat states 

The simple JC model is not, however, always valid. For 
example, if the atom's kinetic energy is of the same order 
of magnitude as the atom-cavity interaction energy, the 
dynamics is significantly changed |lCt ITU . Thus, for very 
cold atoms, the kinetic energy term for the atomic center- 
of-mass motion must be treated quantum mechanically. 
In the standard JC model, the atom is either assumed 
to stay still relative to the cavity mode, or to have a 
large kinetic energy; in both cases the atomic motion is 
described classically and the kinetic-energy term may be 
excluded from the Hamiltonian. Another simplification 
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of the JC model is that the spatial shape of the cav- 
ity mode is not taken into account and it is assumed 
constant. This is, of course, not always valid, since an 
atom traversing a cavity will see a mode that varies with 
respect to the atomic position. The mode variation is 
given by the particular shape of the electric field and 
is, therefore, space dependent. The proper approach in 
such a case is to introduce an atom-field coupling that 
is position-dependent g{x) — d ■ E{x)/Ti, where d is the 
dipole moment of the atomic transition and E{x) is the 
electric field of the cavity mode involved. For a smooth 
coupling and small velocities, the atoms see the cavity 
mode as an effective potential which is, in the adiabatic 
limit, given by ±-\/A 2 /4 + g 2 (x), here A is the atom- 
cavity detuning. Consequently, the atom experiences an 
effective force from the potential and it ma y, fo r instance, 
be reflected or transmitted 
even trapped inside it 
situation in which the atom experiences both an effective 
cavity potential and an external potential has also been 
discussed [22j. Today it is possible to trap ions inside 
cavities even using external traps |23| . which open up 
new possibilities for realizing certain desirable systems. 

The shape of the cavity mode depends on the bound- 
aries of the cavity, the most commonly considered being 
Gaussian, standing wave, travelling waves in ring cav- 
ities, and whispering-gallery modes. We consider here 
a standing wave mode with a wave number q, g(x) = 
g cos(qx), where g is the scaled strength of the cou- 
pling. For such a system, an extended JC model includ- 
ing atomic centre of mass motion and a standing wave 
coupling has been studied in a large number of papers, 
here we just mention a few. Th e dy namics has been anal- 
ysed in, for exainplejpilllilill^ 
while in |30t I3TL I32L l33| . approximation methods are used, 
such as Raman-Nath, Bragg, tight-binding or large de- 
tuning. In the Raman-Nath approximation, the kinetic 
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energy term is neglected, and this has been assumed in 
Refs. 0) HH Hil- where effects from various measure- 
ments on the field or the atom have been studied. 

Clearly, with a standing wave coupling, the Hamilto- 
nian is periodic with period A = 2?:/q. The spectra of pe- 
riodic Hamiltonians are known to consist of allowed ener- 
gies in forms of bands, separated by forbidden gaps. They 
are most commonly treated using the Floquet theory, 
which has been done, for example, in Refs. pil l31j . 
An interesting observation is that the Brilluin zone is now 
twice as wide as in the usual case for one dimensional pe- 
riodic Hamiltonians. This derives from the fact that the 
two atomic levels are coupled to the motion, in contrast 
to electrons in solid crystals where electronic spin flips are 
not coupled to the lattice potential. In an atom-cavity 
system, every time the atom absorbs or emits a photon 
and gets a momentum 'kick', its internal state |±) is also 
changed. Hence, in the rotating wave approximation, an 
emission (absorption) must take place between two ab- 
sorptions (emissions). The symmetry of the system is, 
therefore, generated by displacement of half the spatial 
period accompanied by an atomic inversion (performing 
this twice yields the spatial periodicity of the Hamilto- 
nian), which renders the Brillouin zone (in momentum 
space) twice as wide. 

In solid state physics, it has proved useful to describe 
an electron propagating within a periodic structure in 
terms of a dispersion relation E = E u (k) where k is called 
the quasi-momentum of the state and v is an index for 
the electronic band; in this picture the electron is consid- 
ered to move freely, with its propagation characteristics 
defined by the dispersion relations. Here, as below, we 
set Planck's constant Ti = 1. If the electronic wave func- 
tion is represented by a (Gaussian) wave packet centred 
around (quasi-) momentum k$, its propagation velocity, 
i.e., group velocity, is given by v g = dE v (k)/dk\k=k - 
Thus, the velocity of the electron is no longer deter- 
mined by its original mass, such that v = k/m, but by a 
renormalized mass defined as mi = k/v s , which depends 
on the dominant quasi-momentum of the state. As the 
tangent of the dispersion curve gives the group velocity, 
which is related to the free space velocity through the 
ration m/mi, the curvature determines the amount of 
spreading of the Gaussian wave packet and likewise it 

defines another mass parameter m-i — ■ In this 

paper we study the dynamics of a two-level atom inter- 
acting with a standing wave cavity mode and discuss the 
effect of masses mi and rri2, replacing the original free 
mass of the atom. 

In ordinary QED, the assignment of mass to electrons 
is an essential part of the renormalization program, where 
the mass may be considered shifted by the presence of 
the zero-point energy of the vacuum; the fact that for- 
mally infinite entities are manipulated does not invalidate 
the general picture. Likewise, one may expect that the 
presence of a finite energy in the field may give its own 
contribution to the renormalized quantities, in particular 



the mass. This quantity is usually considered to be too 
small to have any observable consequences. In a cavity, 
on the other hand, the coupling of an atom to the cavity 
modes is enhanced, and it may be possible to interrogate 
the effect of the field on the mass. 

In most setups for atom-cavity QED experiments, the 
atom is prepared in some initial state outside the cavity 
and is then allowed to propagate through the cavity field. 
Provided that the photonic wavelength A is small com- 
pared with the cavity length, the system can be treated 
approximately as periodic, and the results of the Flo- 
quet theory are applicable. If, for instance, the atom is 
prepared with a kinetic energy that lies in a forbidden en- 
ergy gap, it cannot enter the cavity but must be reflected 
from it, possibly with a flipped internal atom- field state 
|±), as will be shown below. When the kinetic energy 
falls within the allowed energy bands, the atom will tra- 
verse the cavity with a (group) velocity v g = fc/mi. We 
also simulate wave packet propagation in these situations 
using the split operator method. The results obtained 
from the Floquet theory and the wave packet propaga- 
tions are compared. Since the mass parameters mi and 
m2 depend on the effective coupling go^/n, where n is the 
photon number, a measurement of mi or mi also yields 
indirectly the photon number inside the cavity. The re- 
flection and transmission of atoms against the cavity may 
also be used for state preparation or 'Stern-Gerlach' type 
of measurements between different internal orthogonal 
states. 

The paper is outlined as follows: First, in Section ITT1 
the Hamiltonian describing the dynamics is introduced 
and solved numerically for the eigenenergies and eigen- 
states in accordance with the Floquet theory. The bare 
and dressed states of the system are presented and an- 
alytic approximations for the lowest band is obtained. 
The more illustrative approach of wave packet simula- 
tions is considered in lllll and the effective masses mi and 
ni2 are defined. Both the propagation of Gaussian bare 
and Gaussian dressed states are discussed and, in order 
to get a deeper understanding of these two cases, we an- 
alyze the comparison between bare and dressed states. 
The masses mi and 7712 are extracted numerically and 
compared between bare and dressed states wave packet 
propagation and also with the masses obtained from the 
Floquet theory. Further it is shown with simulations how 
atoms may be reflected or transmitted by the cavity mode 
and we discuss possible applications for state preparation 
and state measurements. Finally, in II VI we conclude with 
a discussion of the results and possibilities to observe the 
mass in realistic experiments. 



II. THE FLOQUET APPROACH 

We describe the atom-cavity system with a Jaynes- 
Cummings model Q that takes into account two atomic 
levels, coupled to a single field mode in the rotating wave 
and dipole approximations; two essential parameters are 
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the atom-field coupling go and the detuning A between 
the atomic transition frequency u> and mode frequency 
ft. Moreover, the field mode is assumed to be a stand- 
ing wave along the cavity axis x, and the parallel atomic 
motion is quantized. The spectrum of the Hamiltonian 
is obtained using the Floquet theory, and it has a band 
structure with Brillouin zones twice as wide as for one- 
level particals |24l 1271 1.3 ll ] . The effects of the band and 
band-gaps will be discussed in Section lTTTl where we con- 
sider physically realistic situations. 



A. Jaynes-Cummings model for a moving atom 

The Jaynes-Cummings model Q describes the interac- 
tion between a two-level atom and a single field mode. As 
mentioned in the Introduction above, the atomic center- 
of-mass motion has been ignored and the parameters are 
assumed to be independent of the atomic position in the 
original Jaynes-Cummings model; these conditions are 
not, however, always satisfied in realistic atom-cavity ex- 
periment. As the atom traverses the cavity, the shape of 
the coupling will be governed by the cavity mode struc- 
ture. For a standing wave mode, with wave number 
q = 2n/X, the mode is a given by g(x) — gocos(qx). 
In most studies and experiments, the atomic velocity is 
large enough that it can be described classically; thus 
its energy merely adds a c-number to the Hamiltonian. 
In such situations, assuming the atom to be point-like, 
the position operator x can be replaced by the classical 
center-of-mass pos ition moving with the velocity v, i.e., 
x = vt, see pffll l4f|. For cold atoms, however, the center 
of mass motion must be considered quantum mechani- 
cally EHm and the kinetic energy operator term must 
be included in the original Hamiltonian. In many of the 
studies where the kinetic energy term is included, how- 
ever, the system is simplified by adiabatic elimination of 
the excited state in the limit of large detuning [4 ll |42| . 

The extended Jaynes-Cummings Hamiltonian, with 
standing wave mode structure and quantized atomic mo- 
tion, becomes 
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Here the tilde notation (~) indicates the original dimen- 
sional variables. Also m is the atomic 'free' mass, capital 
P and X are the momentum and position operators, a 
and boson annihilation and creation operators for the 
cavity mode, and the a-operators are the Pauli matri- 
ces acting on the internal two states of the atom. Note 
that we only consider quantized motion in one dimen- 
sion along the cavity axis. The coupling will be taken as 
g(X) — 2g~QCOs(qX). The evolution Hamiltonian is writ- 
ten in the interaction picture with respect to the 'free 
Hamiltonian' Ho ~ hQ(^as + a^a + |) and A = uj — Q is 
the atom-cavity detuning. 

Since the Hamiltonian has been given in the rotating- 
wave approximation, the total number of excitations in 



the system N = a^a + ^03 is a conserved quantity and 
the dynamics therefore splits up into separated decoupled 
subsystems for each number of excitations. In the joint 
Hilbert space of the internal atomic state and the cavity 
mode, we define the basis states as 
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= I l,n), 



where the atomic states | j) and | |) are the atomic upper 
state and lower states of the transition, and \n) are the 
cavity mode Fock states. Using this basis and scaled 
parameters, the Hamiltonian assumes the form 
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(3) 



where, after introducing a characteristic time and length 
scale T s and X s , the scaled parameters are expressed in 
terms of the old ones according to 



9 = T s g, 



T s A = T s (u-(l), 



(4) 



and n is the photon number. We will take the photon 
momentum q to define the characteristic length scale as 
X s = 1/q. We shall consistently indicate q in all equa- 
tions below, but use the numerical value q = 1 in all the 
figures in accordance with the chosen length scale. In 
this way, momenta k will be given in units of q and the 
relevant parameters of the model are go/q 2 and A/q 2 . In 
most of the following analysis, we will assume the atom 
to be initially in its ground state | J.) and the mode to 
contain one single photon. We should emphasize that the 
Hamiltonian then becomes identical to the on describ- 
ing a two-level atom interacting with a classical standing 
wave field. Therefor, our model may not only be used 
for describing atom-cavity QED dynamics, but also the 
interaction between two-level atoms and classical fields, 
for example, if the cavity is driven by a classical source 
or the field is given by a laser beam. More general sit- 
uations of the quantized field could be considered in a 
straightforward generalization, but in this paper we only 
discuss the basic features and keep the model as simple 
as possible. 

An interesting observation is that for zero detuning 
A = 0, the unitary operator 



n/2 



1 1 

1 -1 



(5) 



decouples the system into two ordinary one-dimensional 
Schrodinger equations with potentials V± (x) = 
±2cos[oa;); these equations are known as Mathieu equa- 
tions 431. 
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Due to the spatial periodicity A of the cavity mode, 
the operator 



T 



(6) 



with A = 2n/q, commutes with the Hamiltonian J2J; this 
symmetry property is the background for the Floquet 
theory. Another, slightly less obvious symmetry is asso- 
ciated with the operator 



(7) 



that also commutes with the Hamiltonian |2J|. This is 
a 'half-period' displacement combined with an atomic 
inversion and it includes the first symmetry operation 
since T — I 2 . Consequently, the first Brillouin zone is 
within —q<k<q, and not within ±q/2, as implied by 
T alone. Physically this derives from the fact that ev- 
ery absorption or emission of a photon flips the internal 
state |±) — * |=f), while after absorption + emission (or 
vice versa) the internal atomic state remains unchanged. 
The center-of-mass momentum in the two step process 
will either be the same or shifted by ±2g, depending on 
the direction of the emitted/absorbed photons. Note that 
the first Brillouin zone has occasionally been defined to 
extend within — | < k < |, which produces two sets 
of dis per sion curves, one for each internal state |±), see 
Refs.MEl. 



Each dressed state is assigned to some energy band (Bril- 
louin zone) v, which is numbered 1, 2, 3, ... for increasing 
energy, and is also indexed with a continuous variable 
fc, which is now called the quasi-momentum; the entire 
quantum states contain all momenta k + Liq. The func- 
tional dependence of the energy eigenvalue on the quasi- 
momentum E u '(k) is called the dispersion curve, assigned 
to each Brillouin zone. The dressed states for each quasi- 
momentum are obtained by solving the secular equation 
given by the infinite matrix Hamiltonian 
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The eigenvalues and eigenvectors of this infinite Hamilto- 
nian are not known in the general case but approximate 
analytical results may be found, see Refs. [30L l3ll l32l l33| . 
For example, one interesting approximation, related to 
the Raman-Nath limit, is when the g 2 -terms are ne- 
glected and it may be solved analytically. Here we will 
not discuss these approximations, but first solve the prob- 
lem numerically and then make different perturbative ex- 
pansions for the eigenvalues. 



B. Energy band structure 

Due to the conservation of the total excitation number 
N = a^a + i(T3, the system Hilbert space H may be 
represented as a direct sum of Ti.i and TL2 that are not 
coupled by the Hamiltonia. We limit the discussion to 
TLi that is spanned by the bare states 



\k + fiq)\- 
\k + m )\- 



\i even 
\x odd 



(8) 



which are energy eigenstates of the system in the ab- 
sence of interaction, with their energies given by £^ — 
2^-(fc + /iq) 2 — (— l) Mi j- The quasimomentum k is here 
limited into the first Brillouin zone — q < k < q and 
the integer index /i denotes the Brillouin zone or, equiv- 
alently, the energy band. Thus the physical momenta of 
the bare states have well-defined values k + \iq and, in 
particular, the internal |— ) state with zero-momentum is 
given by \ipo(0)) while the |+) state, by |t/>_i (<?)). 

The energy eigenstates of the interacting Hamiltonian 
given by Eq. 



H\<j> v (k)) = E v (k)\cf) u (k)) , v = 1,2,3,. 



(9) 



are called dressed states and they can be expressed as 
linear combinations of the bare states 



\M*>))= E <(fc)|^(fc)>. 



(10) 




fi— — 00 



FIG. 1: The lowest lying bands of the Hamiltonian given by 
Eq. @ for the first Brillouin zone. The dressed energy bands 
are marked on the y-axis with the dominant bare state index 
(i. Crosses shows the bare energy bands for excited states |+) 
and diamonds bands for ground bare states |— ). In the last 
figure (d), the coupling is so strong that the bare and dressed 
energies starts differ considerable. The parameters are given 
on top of each figure. 

In order to solve the problem numerically, the Hamil- 
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tonian has to be truncated at some dimension n. For 
small k and relatively low bands, this n should be cho- 
sen odd, in order to be consistent with coupling to equal 
number of states in both 'directions' from a given initial 
state. For n = 1, we obtain the bare eigencnergy, for 
n = 3, the bare state |^g) couples to the states 
and so on. In Figs. ^ ( a ) _ (d) we show the lowest-lying 
bands for the first Brillouin zone, obtained numerically 
with n = 201, for the parameters (a) A = 0, (b) A = 1 
and (c) and (d) A = — 1. In (d) the coupling is 20 times 
as large, making the bare and dressed energies differ sig- 
nificantly. On the y-axis the dressed bands are labelled 
with the corresponding dominant /i-value that the bare 
eigenenergies would have had in the limit of weak cou- 
pling, and diamonds indicate a bare lower state |— ) en- 
ergies and crosses bare states |+) energies. When adding 
the coupling, the crossings become 'avoided'. Note how 
the gap size decreases with the band index v 1 indicating 
that the state is more weakly coupled to far-away lying 
states. The crossings between even-even or odd-odd 
/i are called Bragg resonances and between odd-even or 
even-odd fi Doppleron resonances |27j . 

In Fig. |21 we illustrate how the presence of the peri- 
odic interaction couples the momentum eigenstates (bare 
states) into dressed states. In Fig. |2Ia), the coefficients 
for the first four dressed states [y — 1,2,3,4) are plotted 
as a function of n for zero quasi-momentum k = 0. Note 
that the solutions are either odd or even in [i and that 
only the first one is not 'degenerate' since all other are 
centered around a crossing. Figure |3fb) shows the same 
coefficients for a non-zero quasi-momentum fco = <z/4, 
and the solutions are no longer symmetric around /i = 0. 
Note that, if any of the coefficients c 1 ^ (for each v) in 
Eq. lfTU|) has an absolute value close to unity, the pres- 
ence of the periodic coupling only modifies the properties 
of a bare state without too much coupling to other bare 
states; this is usually the case away from the crossings. 
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FIG. 2: The expansion coefficients for the four lowest 
dressed states v = 1,2,3,4, Eq. 11 PL of bare states /i for 
(a) k = and (b) k = 1/4. Black bars corresponds to even 
fi's with ground state atom, and white bars to odd fi's with 
excited atom. The coefficients are symmetric around fi — 
only for fco = 0. The parameters are in both plots the same 
as for Fig. □ A = and g = 0.05. 



C. Extraction of the effective parameters 



fraction-like expression 



We now look for an analytical expression for the energy 
eigenvalues of the Hamiltonian (|11|) and for the dispersion 
curves. Denoting the bare-state energies by 



(k 



2 1 ' 2' 



(12) 



the energy eigenvalue equation H\(f>) — E\4>) (for a pre- 
defined k) can be written as a recursion equation 



= Ec„ 



(13) 



This equation has, naturally, an infinite number of solu- 
tions corresponding to different bands. Truncation of the 
recursion symmetrically around some index v and elimi- 
nation of the expansion coefficients yields a continued 



E = S V 



which has a form of an iteration equation. We assume 
that this equation is most applicable in the region where 
one bare state dominates the dressed state and the trun- 
cation of the Hamiltonian (|ll|l or, equivalently, of the re- 
cursion equation (|13|l is performed symmetrically around 
this state. Note that the truncation of the Hamiltonian 
is also an effective expansion in the coupling constant 
go since, for a small coupling, each base state only cou- 
ples to a few neighboring bare states while, for a large 
coupling, many bare states are required to represent a 
dressed state. This is not, however, true in the vicin- 
ity of level crossings where two bare state couple to each 
other even over many intermediate states. 
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We illustrate the truncation error with 
S(g ,A;n) = \E{k) - E n (k)\, 



(15) 



where E(k) represents the exact dispersion curve and 
E n (k) the one obtained from a truncated n x n Hamilto- 
nian. Wc emphasize that the error depends only on the 
parameters go and A due to the chosen length scale; in 
physical units they both contain the photon wave number 
since go oc go/q 2 and A oc A/<j 2 . As already mentioned, 
we use q = 1 throughout the paper. 

Figure |31 illustrates <5(<7o>A;n) as function of go and 
A, for k = 0. An almost identical plot of the error is 
obtained for k = q/4, therefore it is omitted. The error 
increases for large couplings, which is easily understood 
since a large coupling means that the initial bare state 
will couple more strongly to other bare states and the 
dimension of the Hamiltonian must be chosen higher. It 
is also seen that the increasing detuning A, which makes 
the diagonal elements in the Hamiltonian larger, yields 
smaller errors. 




o -S 



FIG. 3: Error estimate <5(<?o,A;n) for n = 5, as defined in 
Eq. 11511 . The quasi-momentum k = 0. The errors for the 
k = 1/4 case is almost identical. 

The number of ^Q-terms included in the continued frac- 
tion (|14|) is related to the truncated size of the Hamilto- 
nian as (n — l)/2. Taking any initial value of E and 
iterating Eq. (jT4"|) it is expected that the iteration con- 
verges to some eigenvalue E v (k) close to the initial value 
of E. For example, for moderate couplings and away 
from crossings, starting with E coinciding with a bare 
energy eigenvalue, the iteration is supposed to converge 
to the corresponding dressed energy eigenenvalue. Thus, 
analytical approximate results are obtainable by trunca- 
tion the continued fraction to some (n — l)/2 terms and 
iterate it j times. From numerical investigations of the 
validity of these two approximations, we draw the conclu- 
sion that for certain rangers of the parameters, especially 
for small couplings go and large detunings A > 1, the 
order of approximation does not need to be very high 
away from crossings. It has turned out that truncating 
the Hamiltonian to a 5 x 5-matrix and iterating the con- 
tinued fraction (|14f> twice gives an eigenenergy, which is 



valid for a large range of parameters. Having an ana- 
lytical expression for the energy, we can easily calculate 
the mass parameters. Below we will give only the ana- 
lytic expression for the case when k w 0, but the same 
procedure could be carried out for other cases as well. 
Since we have assumed |fc| <C q, we expand the eigenen- 
ergy around k = 0, but we will also expand the result in 
powers of <?q. The result obtained becomes is given by 



A 

' 2 



,4gg , 

<j 2 +2A 

16q 2 ffg 
(<J 2 +2A) :i 



4(7g 2 -2A)gg 
^(g a +2A) 



O(g e o) 



4(lll<f-46Aq 4 -28A 2 q 2 -8A 3 ) fl 



9 4 (g 2 +2A) 5 



+O(g 6 o))k2 + O(k i ). 



(16) 



III. PROPAGATION OF WAVE PACKETS 

In this Section, we analyze the propagation of Gaussian 
dressed and Gaussian bare states using the effective mass 
parameters and compare the results with wave function 
simulations. We also discuss the physical difference of 
initial bare and dressed states. 

Propagation of an initial Gaussian state can be un- 
derstood in terms of effective parameters, such as the 
group velocity and the effective mass, which depend on 
the dispersion curve and are evaluated at the dominant 
quasi-momentum fco of the wave packet. We assume that 
the wave packet inside the cavity is initially described by 
a Gaussian momentum wave function 



<p(Js) 



1 



(17) 



where A^ is the width of the momentum distribution and 
it is related to the initial width of the position distribu- 
tion according to A& = 1/(2?: A x ), whence the initial 
state is a minimum uncertainty state. Within its range, 
the energy of the dressed states (belonging to the band 
v) can be expanded as 

E(k) « E(k ) + v s (k Q )(k - fc ) + (A - fco) 2 

= E ° + k ° + ^TTM fc o(fc - fco) 

+ 2^j( fc -M 2 , 

(18) 

with Eq = E(k = 0). We have chosen to use mass pa- 
rameters for each term in the expansion. Note that, for 
a free particle, this is merely an expansion of the energy 
term E = ^fc 2 around fco, and all the mass parameters 
vci\ coincide with the natural mass of the particle. 

The three mass parameters in Eq. (|18H can now be 
given the following interpretations: mo yields the en- 
ergy of a dressed eigenstate as -E'(fco) — E = 2m^^o 
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and, hence, it defines the phase velocity of the dressed 
state as v p — E Q /k Q = 2^-fco> m i defines the group- 
velocity-quasi- momentum relation as fco = miv s , and 
mi determines the mass associated with the wave-packet 
spreading. Note that, even though we do not consider 
the case here, TO2 coincides with the conventional effec- 
tive mass m* that determines the acceleration caused by 
an external force acting on the particle wave packet. The 
first mass parameter mo is, however, nonphysical since its 
value depends on the choice of zero energy level. It may, 
however, affect some interference experiment, but we do 
not consider it here, since it is not expected to effect the 
propagation. Equation l|16() can be used to deduce ap- 
proximate values for mi and TO2 near zero-momentum. 



Gaussian throughout the evolution, while the initial bare 
state splits up in three main sub-packets corresponding 
to the bare states \ — ,k ), \ + ,k — q) and | + ,fc + q). 
In order to prepare the initial dressed states, the cou- 
pling amplitude is taken to be time dependent <?o(^) and 
is turned on adiabatically from go = up to the final 
value go(t) = go. In that way, the state remains dressed 
during the turn-on assuming an adiabatic switch on. As 
go(t) has reached the final value go, the wave packet has 
already broadened, so after the preparation process we 
do no longer have A^, = l/(2Afc). The numerical method 
used for wave packet simulations will be discussed in Sec- 
tion nun] 



A. Propagation of Gaussian dressed states 

By Gaussian dressed states we mean wave packets cen- 
tered around the quasi-momentum fco, where each wave 
component belongs to the same energy band of the inter- 
acting Hamiltonian, 

/oo 
^)|<Mfc))e- i£ " Wt dfc. (19) 
-oo 

In principle, the integral should be limited to the first 
Brillouin zone (or to any one Brillouin zone); in this ex- 
pression it is assumed that fco is sufficiently far from its 
boundaries so that the momentum distribution ip(k) is 
negligible outside the zone. 

With the use of the expansion of the dressed states in 
terms of bare states, Eq. (fTTTf) . and of the integral equality 



JZo^k) e i(»»-Mfc)t) d fc 
(*-v) 2 



(20) 



J(kox-E(ko)t) 



2tt 



(the approximate value derives from neglecting higher 
terms in Eq. I|18|0 . the Gaussian dressed states evolve in 
time as 



\Mko)) 



2tt 



•A (21) 

and has the time-dependent width (t) — | ^ — h | . 
Here we have further assumed that the expansion coeffi- 
cients c^(fc) remain constant within the Gaussian mo- 
mentum distribution. Inclusion of a correction term 
c 1i{k) ~ c jK^o)+G^(fco)(fc— &o) gives rise to an additional 
term in the integral, but it still has a Gaussian envelope 
moving with the same group velocity. In Figs. 0] (a) and 
(b) we show the propagation of initial Gaussian dressed 
and bare states. The dressed state stays approximately 
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FIG. 4: The evolution of (a) an initial Gaussian dressed 
state and (b) an initial Gaussian bare state. The dressed 
state stays localized around an average position, while the 
bare state clearly splits up. The three sub-packets in (b) 
correspond to the bare states | — , fco), | + , fco — ?) and | + , ko+q), 
which is also seen in the inset showing the final momentum 
distribution |(fc|4')| 2 . The parameters are the same in both 
plots, g = 0.001, A = 0, x = 0, fc = 1/4 and A 2 k = 1/10000, 
and in (b) the atom is initially in its lower state. The contour- 
bar shows relative values. 
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B. Comparison between bare and dressed states 

While in most experimental situations an atom enters a 
cavity in a well-defined bare state, the mass parameters 
defined above characterize a dressed state and are not 
directly applicable. We introduce therefore the overlap 
(often called Fidelity) between bare and dressed states 

F^(fc) = |(^|^)| 2 , (22) 

which is the same as the absolute square of the coeffi- 
cients in Eq. (|1U|I For a relative small coupling and 
away from crossings, it is clear which bare state 
has the largest overlap with a given dressed state \<j>„)- 
However, near crossings two bare states will become im- 
portant, but still it is easy to decide which ones. Note 
that, given a quasi-momentum k, the crossings depend 
on the parameter A, and the fidelity will change drasti- 
cally around some particular value on A, see Fig. ^ For 
example, picking k = q/A and A = 0, the largest overlap 
with the dressed state with v = 1 is assumed to be with 
the g. = bare state, assuming a small coupling go, while 
for — 3q 2 /A < A < — q 2 /A the largest overlap is expected 
for the fj, = — 1 bare state. 

In Fig. [SI we display the coefficients mjirn-i) and u g = 
ko/mi as function of go and A. In Figs. (a) and (b) 
1/(7712) is shown for k = and k$ = l/A respectively, and 
in Fig.EJc) the group velocity is plotted for k = 1/4. 

From Fig. it is clear that 1/(7712) differs considerably 
from the unperturbed value m = 1. It is now of interest 
to see how well the corresponding bare and dressed states 
overlap for the same parameter rangers. We calculate the 
overlaps F v=1 ^=°{k = 0) and F y=1 ^=°{k = q/A) for the 
two examples. The results are plotted in Figs. Eta) and 
(b). As justified earlier, a small coupling increases the 
overlap. Note that there are ranges where the masses mi 
(or v s ) and 1112 are shifted significantly from the free mass 
m and still the overlap is large, meaning that an initial 
Gaussian bare state should evolve approximately freely 
with mass parameters mi and mi. In the wave packet 
simulations below it will be shown that the shifted masses 
can be observed also for propagating bare states. 



C. State propagation simulations 

Independently of the chosen basis (bare or dressed 
states), the time evolution of the atom-cavity system 
is governed by a time-independent Hamiltonian H = 
p 2 /(2m) -I- V(x); with an initial state |\P(f = 0)) it is 
written as 

|*(t)) =c- iHt \y{t = Q)). (23) 

For an x-dependent potential V(x), we find, in general, 
that the kinetic energy term of the Hamiltonian does not 
commute with the potential. The exponential can there- 
fore not be separated into one momentum exponential 




FIG. 5: The scale mass parameters 1/(7712) and v g = ko/m\ 
as functions of go and A. The first two plots (a) and (b) 
show 1/(7712) for k — and k — q/A, respectively, and (c) 
the group velocity v s for k = q/A, where the unperturbed 
velocity (go = 0) should be vo = 1/4. The result derives from 
numerical diagonalization of the matrix li lt and numerical 
calculation of the coefficients of the dispersion curve around 
the point k of interest. The bare energies of states g = 0, —1 
cross for A = kq — q 2 /2 and the dominant bare states of the 
lowest band changes. This leads to a drastic change in the 
effective parameters, as seen in the figures. 

and one spatial exponential. However, if the time of 
propagation St is chosen sufficiently small, the error of 
splitting up the exponential becomes negligible 

|*(* + St)) « e- l ^ st e- iV ^ 5t \^(t)). (24) 

This is called the split operator method, see Ref. 1 16] , and 
the procedure goes as follows: Starting with some state 
1^(0)) in the x- representation, we multiply it by the ex- 
ponential exp(— W(x)St), we then take the Fourier trans- 
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the internal state 



the wave function is 




-5 3 




-5 3 



FIG. 6: The variation of the overlap (a) F 1,0 (0) and (b) 
F 1,0 (g/4), as defined in Eq. t\2'3 . with respect to go and A = 
0. For positive k, a huge change in the overlap occurs when 
A = kq — q 2 /2, which was also seen in the previous figure of 
the mass parameters. For such A, the corresponding quasi- 
momentum k is located at a gap, and as A changes more, the 
lowest dressed states v = 1 will have a larger overlap with bare 
states with /i = ±1 rather than /i = 0. This can be clearly 
understood from Fig. 0(b) and (c). The photon momentum 
q, which was used to define the characteristic length scale, is, 
as earlier 1. 



form of the state and obtain it in the p-representation, 
and then we multiply it by exp [— ip 2 8t/{2m)\ , and finally 
transform it back to the x-representation with the inverse 
Fourier transform, which gives us \^>{5t)). Repeating this 
N times we obtain for t = N5t. 

For the system considered here, the kinetic and poten- 
tial terms are 



Ekin 



V(x) 



2 



1 
1 



(25) 



If the atom is initially in a bare state corresponding to 



l*(* = 0)> = 



Jdkip{k)\k) 



(26) 



while, for an initial dressed state, it is given by Ea. H19|) 
for t = 0. 



D. Propagation within the cavity 

For a given kg, the properties of the dispersion curve, 
and hence of the mass parameters, depend on the phys- 
ical parameters of the system. From Eq. (|16fl . it is ex- 
pected that effects are largest for A = as long as fco 
is far from a crossing, and thus we choose A = in the 
analysis below. Furthermore, we choose to calculate the 
quantities for fco = and fco — <?/4 in the lowest band 
v = \. 

We will, in this Subsection, show numerical simula- 
tions and calculations of the effective parameters for a 
Gaussian wave packet inside a cavity, but first we inves- 
tigate the evolution of a bare Gaussian wave packet that 
is prepared inside the cavity with a momentum fco = q/2, 
corresponding to a gap. At an energy gap, the dispersion 
curve has a vanishing first order derivative. Hence, if the 
wave packet has a large overlap with the corresponding 
dressed states, its group velocity is close to zero. In the 
weak coupling limit, the time dependent bare state can 
be expressed in terms of dressed states approximately as 



(27) 



+ e- iB " = *|0 2 (fc o = q/2))) , 
and using that the gap size is approximately 2go it follows 



\W, fco = q/2)) 



l 

V2 



e i ^\ ( f> 1 {ko = q/2)) 
-totfafa = q /2))) , 



-i£"t 



.(28) 

Thus, the bare state wave packet is supposed to oscillate 
between upper and lower states |±). Figure shows the 
evolution of such a prepared state. The upper figure 
shows the wave packets in the x-representation for upper, 
lower and combined internal states, while the lower plot 
displays the momentum packets. 

When the atom propagates inside an infinite cavity, it 
should be possible to describe the atomic wave packet as 
a freely evolving particle with effective mass parameters 
defined in Eq. (|18|l . which will be analyzed next. The 
time of flight t{ between xq and X{ is given by 



X[ — Xq mi (X[ — Xq) 



m initial 



(29) 
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FIG. 7: The evolution of an initial Gaussian bare wave packet 
prepared inside the cavity with A^ = 500 and ko = 1/2, cor- 
responding to a quasi-momentum of the first forbidden energy 
gap. The upper plot shows the wave packets of the ground, 
excited and total (ground plus excited) atomic state in the 
^-representation. The lower figure displays the wave packet 
of the ground and excited atomic state in the momentum k- 
representation. It is clear that we have a spreading in x and 
also a small drift of the lower atomic wave packet. The pop- 
ulation 'Rabi' oscillates between the two states with a period 
Tr si 7r /go, corresponding to the 2n divided by the energy gap 
size 2go in first order. The remaining parameters are A = 
and go = 0.05. The contour-bars gives the relative values. 



where ^initial is the velocity of the atom in absence of the 
cavity interaction. We will simulate propagation of both 
bare and dressed wave packets. In the simulation of bare 
state propagation we use averages, for the atomic state of 
interest, that is for the wave packet of the ground state, 
(x) = ^jzjzt- > in other words we select the wave packet of 



the ground state 



the averages. Physically that assumes that a projective 
measurement is carried out on the lower state before that 
position measurement is done. For the dressed states the 
averages are calculated for the whole atomic wave packet 
and no projective measurement is assumed. We propa- 
gate the initial wave packet from a certain (xo) for a time 
tf and calculate the final average position (xf), and de- 
duce the propagation velocity v — i^lzJ2°l ( which should 
coincide with the group velocity v g = kg /mi. We use 
the parameters of one of the previous examples, A = 0, 
A 2 , = 2500 and fc = <?/4, but vary the coupling strength 
go', the final time is taken to be tf = 1000. During the adi- 
abatic preparation, the dressed wave packet will already 
broaden, whence its initial width is no longer A 2 , = 2500, 
but this is not expected to affect the group velocity. In 
Fig. [S] we compare propagation velocities for dressed and 
bare states with the group velocity obtained from the Flo- 
quet theory. The three curves coincide rather well, but, 
for the case of initial bare state, the agreement is good 
only for the lower initial state |— ), as described above. 
The whole atomic wave packet, containing both upper 
and lower atomic states, gives a value (xf) different from 
the one using just one atomic state. Thus, we empha- 
size again that, extracting the propagation velocity with 
the bare states assumes that the position measurement 
is carried out for the atom in its lower state, in other 
words, some of the outcomes will fail, namely when the 
atom is measured while being excited. The slight differ- 
ence between dressed state propagation and the Floquet 
theory might be due to non-adiabatic generation of the 
initial dressed state. 

This scheme clearly gives a good measure of the group 
velocity and from that it is possible to calculate the mass 
mi. 

The mass parameter m?, enters the wave packet dy- 
namics through the broadening of the Gaussian dressed 
state according to 



Kit) 



i 



Afc 
m 2 



(30) 



and normalize it and then calculate 



Initial Gaussian bare states do not, however, remain 
Guassian but split into several wave packets, correspond- 
ing to a number of dressed states. In Fig.|H(b) the initial 
bare state splits into three main wave packets and, even 
though the two outer ones may be small, they contribute 
significantly to the variance since they are far from the 
center packet. Thus, for the calculation of the bare state 
variance, the packets in the tails |x| > 150, which mostly 
corresponds to excited atomic states, are excluded when 
(x 2 ) is calculated. In Fig. the mass ui2 is displayed as 
function of go. The time of propagation is tf = 10, and 
the mass is obtained by a least square fit of Eq. (fTTTTfl . 
knowing A 2 (t) and t. Note that when we prepare the 
Gaussian dressed states by an adiabatic turn on of the 
coupling amplitude, Eq. (pTTTf> is valid, but the mass m-i 
is not constant during the turn-on, which must be taken 
into account. This is carried out by assuming a constant 
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FIG. 8: The group velocity n g as function of the cou- 
pling, obtained, both from the wave packet simulation, circles 
shows the result from bare state propagation and crosses from 
dressed state propagation, and the Floquet theory (solid line), 
when fco = g/4. A clear agreement is seen, but unexpectedly 
the bare state result coincide better with the Floquet one, 
which might be due to a not fully adiabatic preparation of 
the dressed state. The parameters are if = 1000, A = and 
Al = 2500. 



FIG. 9: The effective mass rri2 as function of the coupling, 
obtained, both from the wave packet simulation of bare (cir- 
cles) and dressed states (crosses) and the Floquet theory (solid 
line), when k = 0. The main reason for the deviation of the 
result of the bare state propagation is that the wave packet 
of the lower atomic state, which is used in the calculation, 
does not stay Gaussian and the tails of the distribution are 
cut off for \x\ > 15 in the calculation of AjJ. The parameters 
are if = 1000, A = and A 2 X = 1500. 



coupling amplitude, adjusted by introducing an effective 
turn-on time t e s. Thus Eq. (|3L)[) becomes 



A£(i) = 



4A^ 



^] (t eS 
m 2 ' 



tf 



(31) 



where the mass is now constant during the turn-on 
process. In lowest order, the mass m-i behaves as oc g^, 
according to Eq. (JTBJl, which is indicated in Fig. and 
it has been confirmed by a polynomial fit to the data, 
where the dominant coefficient is the one of g^. 

The measurements of the masses proposed in this Sub- 
section may be used for state preparation. The masses 
are functions of go, A and q for the one excitation case, 
but in general it is a function of the effective coupling 
go\/n, where n is the photon number, instead of go. Thus, 
for a general initial state of the cavity field, a perfect 
mass measurement will give ?ni,2 = mi^igos/n, A, q) for 
any n. Only masses mijs(goy/ri, A, q) with n's that has a 
non-vanishing probability from the initial photon distri- 
bution, can be detected. Knowing mi or 7712 and go, A 
and q it is possible to solve for n, meaning that measuring 
the mass reduces the field to the Fock state \n), a sort of 
projective measurement. Physically this means that an 
initial state |$(t = 0)) splits up in a set of states |$(i)) n 
for the various photon numbers n. If all overlaps between 
different final states |$(i)) n vanish it is possible to sep- 
arate the individual parts with one single measurement. 
However, having non-overlapping final states seems un- 
likely in realistic experiments. 



E. Scattering and transmission 

In order to simulate atomic scattering and transmission 
by the cavity, the mode function coupling g{x) needs to 
be modified by multiplication of an envelope function 
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tanh 



tanh 



(32) 



This function goes to zero for large \x\ and is centered 
around x = 0, the cavity length is given by 1x\, and 
the 'slope' how fast it turns on/off is determined by x e . 
Thus, the coupling will approximately be zero for x < 
—x\ and x > x\, and in the interval [— afi,a;J it behaves 
as g(x) = 2gocos(qx). We choose x\ such that 1x\ <C 
A = 2n/q, and x e such that the coupling is turned on/off 
smoothly, but fast enough for a non-adiabatic transition. 
The envelope function allows us to test boundary effects 
from not having a completely periodic potential; it turns 
out, however, that the dynamics follows very well the 
Floquet predictions as long as A <C 2x\. 

An incoming atom entering a cavity has been studied 
in a series of papers [l2l IT3I Il4j. were, in particular, the 
transmission and reflection coefficients have been calcu- 
lated. In these papers, either a mesa function or a hy- 
perbolic secant squared coupling g(x) was used to obtain 
analytical solutions. The spatial wave function was as- 
sumed to be a plane wave and not a propagating wave 
packet. In the case of a mesa function or a hyperbolic 
secant, the atom effectively sees a potential barrier or a 
well with an amplitude ±-y/A 2 /4 + g 2 (x)n, where n is 
the photon number. The case of a standing wave shape 
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of the coupling is, however, different; the spectrum of 
allowed energies is then formed by bands, with forbid- 
den gaps in between. Thus, an atom entering the cavity 
mode must have an allowed energy in order to traverse 
the cavity. If the atom has an energy falling into for- 
bidden energy gaps, it must be reflected (assuming the 
tunneling rate through the cavity to be negligible) . 

Going back to Fig. ^ we see that, when A = 0, the 
crossings occur for (fi + l/2)q and /zg, where, as before, 
/i = 0,±l,±2,... . We know that these become avoided 
crossings when the coupling is turned on, and they give 
the forbidden energy gaps in the spectrum. The gap size 
decreases for increasing band index u, and the largest 
gap is at k — ±g/2 between the first two bands. If the 
momentum distribution ip(k) of the incident atom is such 
that (k) » 0.5g and the energy spread, due to the spread 
Afc, is smaller than the gap width, the atom should be 
reflected when it approaches the cavity mode, and from 
energy conservation we must have (k) w 0.5q =>■ (As) ~ 
— O.bq. The atom thus shifts the momentum by one unit 
q and, due to the form of the Hamiltonian (0, also the 
internal atomic-cavity state flips. Hence the quantum 
state is changed, by the reflection against the lowest gap, 
as 

J dk(p{k)\k)\±) =*> J dkip{k)\k~q)\T). (33) 

On the other hand, the internal atomic state will not flip 
for scattering against a gap corresponding to momentum 
fiq. If the atom goes from | f) (| J,)) to | J.) (| |)), the 
photon number must also change as n — > nil. This 
means, for example, that starting with the vacuum and 
reflecting N excited atoms from the cavity allows one to 
create an TV-photon Fock state in the mode. If instead 
the atoms scatter when they are initially in their ground 
state, the mode will be cooled down towards lower photon 
numbers. 

Figure ITUI shows the results of a numerical simulation, 
in which an atom reflects from a cavity. The atom is 
initially in its lower atomic state |— ) and ends up in |+). 
The cavity starts at around x = —150 and ends at x = 
150, x\ — 1500 and x c — 50, and the flip takes place at 
the edge of the cavity. The scale of the contours is chosen 
such that even the small amplitudes of the wave packet 
is seen; a closer look on the scale indicates that almost 
all the population is reflected. 

If the energy spread of the incident atom is broader 
than the band gap, the tails of the packet are in the 
allowed region of energies, and can therefore enter the 
cavity. Figure ITT1 depicts the same system as Fig.^J but 
the momentum distribution is much wider, hence part of 
the wave packet is transmitted. This phenomenon can 
also be observed on the momentum representation of the 
wave packet: The part of the wave packet having al- 
lowed energies is transmitted, while the forbidden part is 
reflected, with a flipped atomic state | [) — > | |). Thus, 
due to the forbidden energies, a 'hole' is seen in the mo- 
mentum distribution. This suggests a way of measuring 
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FIG. 10: Scattering of an atom against a cavity. The cav- 
ity is located within \x\ < x\ for x\ = 1500 (marked with 
dashed line) and x e = 50. In (a) we show the ground state 
atomic wave function |( — | (a; | \I/ ) j 2 and in (b) the upper state 
|(+|(a;|\t)| 2 as function of position x and time t. Note that the 
atom starts to propagate in its lower state at xo — —2500 at 
the time t = with momentum ko = q/2 until it reaches the 
cavity boundary at x = —1500. It is reflected and its internal 
state is flipped to the upper state | ]). The reflection is due to 
the fact that its energy coincides with an energy gap; hence 
the atom absorbs the photon from the cavity field and the 
impact changes the direction of propagation. Note the scale 
of the contours, which indicate that almost everything is re- 
flected, as confirmed in Fig. I13l for the atomic inversion. The 
parameters are here are go = 0.01, A = 0, and A 2 = 2500. 
The contour-bar shows relative values. 

the band/gap structure, since the energies of the reflected 
atoms correspond to gaps. The momentum wave packet 
for the |— ) atomic state is shown in Fig. We observe 
that the momenta close to k = q/2 are 'burned out'. 

The atomic inversion, defined as (era), indicates the 
amount of the atomic population that is in the upper 
respectively the lower state. If (173) = —1, all the popu- 
lation is in the lower state, while if (173) = 1, the atom 
is completely excited. In Fig. 1131 we show the inversion 
for the two examples above, = 2500 and A^ = 100, 
as function of time. Here it is obvious that, in the first 
case, almost the entire population is flipped, while in the 
second most of it is not flipped. We note that, after the 
atomic wave packet has reached the cavity, an oscillating 
behavior can be observed. The reason is the interference 
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FIG. 11: This figure shows the same as Figs.^](a) and (b), 
except that now the momentum distribution is much wider, 
AjJ is chosen 100 instead of 2500 as in the previous plot. The 
consequence is that energy spread exceeds the band gap and 
the tails of the momentum wave packet are transmitted rather 
than reflected. The contour scale is the same as in Fig. fTTTI 
and a closer look shows that a much larger part of the wave 
packet is propagating inside the cavity. Note how the wave 
packets spreads much faster in these plots. 



between the dressed states involved. 

So far the atom has entered the cavity with a momen- 
tum (k) that lies within a band gap, but we now consider 
the situation in which it is within the first allowed band 
v = 1. We still choose the atom to start in its ground 
state |— ), but now with fco — <?/4, halfway between k = 
and the gap at k = O.bq. If the overlap between the ini- 
tial bare state and the dressed state \<j>^) is close to unity, 
the population of the internal states of the system stays 
approximately constant, and the atom will traverse the 
cavity without changing between its two states | f) and 
| I) nor change its momentum state. It will, however, 
experience a shift in the masses m\ and ?7i2- The over- 
lap for these parameters is F"= 1 ^=°{k () = 1/4) = 0.998, 
which is close to 1. The result of the simulation is given 
m Fig. [H The atomic inversion is shown in Fig. 1151 from 
which we see that ~ 99.8 % of the population is in the 
lower state during der the situation in which it is within 
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FIG. 12: The momentum distribution of the lower atomic 
state, |( — |(fc|4')| 2 , corresponding to Fig. 1111 as a function of 
time t. The forbidden energies are 'burned out' when the 
packet hits the cavity. The contour-bar gives relative values. 
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FIG. 13: The atomic inversion (0-3) as function of time t for 
the example in the previous Fig. 1101 (solid line), and Fig. 1111 
(dashed line). It is clear that at the cavity boundary, the 
atom flips from | J.) to |) in the first example, but when 
the momentum width becomes larger than the gap, parts of 
the wave packet will not be reflected and therefore will not 
be flipped. 



the first allowed band the time inside the cavity, as pre- 
dicted by the overlap. After leaving the cavity, the atom 
regains its original internal and momentum state. 

As mentioned above, scattering N excited atoms 
against an initially empty cavity, so that they are all 
reflected in the lower state, creates a Fock state with ex- 
actly N photons. Other interesting situations may be 
realized as well with the same idea. Assuming the cavity 
to be empty from the beginning and the atom in some 
linear combination of upper and lower state and with 
a momentum corresponding to the first gap, the state 
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FIG. 15: The atomic inversion (cr z ), corresponding to Fie. 1141 
Note that when the atom has traversed the cavity, it regains 
its original state. The asymmetry is due to broadening of the 
wave packet while traversing the cavity, leading to a longer 
transit time as it leaves the cavity compared to when it enters. 



a model for creating entanglement between the involved 
degrees of freedom, momentum, atomic and field states. 



FIG. 14: Transmission of the atom through the cavity. The 
plots are the same as in Figs. llOf a - ) and (b), but now the 
initial momentum is ko = q/4, and the atom energy falls 
on the first allowed energy band, contrary to Figs. I1UI and 
1111 The overlap between the initial bare state I^q 0-1 ^ 4 } and 
the dressed state l^ " 174 ) is 0.998. Thus, the propagation 
can be described using the effective parameters and the wave 
packet remains Gaussian; only a small amount of population 
is transferred into the |+)-state and the Gaussian wave packet 
spreads according to the effective mass m,2 and traverse the 
cavity with a velocity v s . The contour scale is the same as in 
Fig. ^] and the cavity boundaries are marked by the dashed 
lines. 



evolves as 

(a\t) + b\ |))k/2)|0) => (a|-«/2,l) + %/2,0))|i). 

(34) 

Since a lower state atom will traverse unaffected through 
the cavity, while the upper state atom is reflected, this is 
realizing an atomic 'Stern-Gerlach' type of measurement. 
If instead, the field is initially in some linear combination 
of vacuum and the one photon state and the atom in the 
lower state it follows that 



i, 9 /2)(o|0) + 6|i; 



(all, 9/2) 



6|T,-?/2>)|0>. 

(35) 

Thus a photon 'Stern-Gerlach' type of measurement is 
realized, assuming that the field is restricted to the states 
|0) and |1). Obviously, the scheme may not be used only 
for separating atomic or field stats, but it works also as 



IV. CONCLUSION 

We have analyzed theoretically the problem of a two- 
level atom propagating within a quantized standing wave 
electromagnetic field, with an emphasis on the mass pa- 
rameters in atomic cavity QED systems. We have ap- 
proached the problem with two different strategies: the 
Floquet theory for the periodic Hamilton operator, and 
the wave packet simulations. The first is a purely math- 
ematical procedure, which offers numerical and approxi- 
mate analytic expressions for the mass parameters, while 
the second is a more physical attack. The analytic re- 
sults offer a hint of the dependence of the masses on the 
physical parameters. Roughly speaking, a small detun- 
ing A, small photon momentum q, and a large atom-field 
coupling go yield large mass shifts. 

The dynamics of an atomic wave packet inside (or en- 
tering) the cavity can be understood in terms of the dis- 
persion curve E v {k) : where k is the quasi-momentum of 
a generalized plane wave and v is the band index. In 
the Taylor expansion around some quasi-momentum kg 
each term defines a mass paramter; in this paper we have 
considered the mass m\ = hko/v g related to the group 
velocity v g of Gaussian wave packets and the mass mi 
associated with the spreading of the wave packet: 



1 

mi 



1 dE v {k) 



h 2 ko 



Ok 



and 

k=k ™2 



1 d 2 E u (k) 



dk 2 



h—k,Q 

(36) 

Note that m 2 is the same as the customary effective mass 
m* introduced to describe, for example, acceleration of 
electrons within a crystal under an external force. These 
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parameters vary essentially over different scales of the 
coupling, as seen in Fig. [SJ Even an arbitrary small cou- 
pling g results in strong effects in the atomic dynamics. 
For g > i, both the group velocity and l/m2 approaches 
zero and finally saturates towards this value, which in 
unsealed quantities can be written as 

Ti 2 a 2 

E- m t =n~g< -7^- = ^kin. (37) 

Hence, whenever the characteristic interaction energy ex- 
ceeds the characteristic kinetic energy of the lowest band, 
the band is essentially flat and the energy becomes nearly 
independent of the quasi-momentum. Using typical val- 
ues A = 1000 nm. and m = 100 a.m.u. we get a char- 
acteristic time scale T s ~ 0.4 [is. giving us the unsealed 
coupling g = g/T s w 10 kHz for saturation. Note that 
dividing our dimcnsionless detuning A and the coupling 
g by T s gives real physical quantities. 

The effective masses are only defined for dressed states 
with well defined momenta, and in many situations, the 
system starts out in a bare state. It is, however, possi- 
ble that the bare states have a large overlap with some 
dressed state and may be assigned an effective mass. This 
has been investigated in great detail in the present paper 
and numerical simulation of propagation for both Gaus- 
sian dressed and Gaussian bare states has been analyzed 
and compared with both each other and with results 
given by the Floquet theory. Both the group velocity and 
the effective mass mi have been extracted numerically in 
the various cases with a good agreement as compared 
with the values obtained from the Floquet theory. 

In most experiments today, the atom and field are pre- 
pared separately and the atom is shot through the cavity 
with a preselected velocity. The state of the complete 
system is then, most likely, not in a dressed state, but 
rather a bare state. We have not, in the text, discussed 
in detail how one may prepare initial dressed states. A 
standard method is to use a slow adiabatic change of 
some external controllable parameter. For the JC model, 
bare and dressed states becomes identical in the two lim- 
its go — > or A — > ±00. By letting the atom traverse a 
Fabry-Perot cavity almost perpendicular to the standing 
wave mode it is possible to have very small initial atomic 
velocities hko / m along the standing wave mode even for 
not so cold atoms by varying the angle of injection; the 
atoms will have a high transverse velocity and a slow ve- 
locity in the ^-direction along the standing wave mode. 
Thus, it is enough to quantize the atomic motion in just 
the x-direction, while the y and z-directions are treated 
classically. In the case of a Fabry-Perot cavity, the trans- 
verse Gaussian mode shapes can then be described by 
the time envelope exp(— t 2 /At 2 ). If the turn-on of the 
transverse Gaussian mode is slow enough, the atom will 
traverse the cavity in a dressed state. Another feasible 
way to prepare the dressed states is to use an external 
ion trap to confine the ion inside an off-resonance cavity, 
and then slowly tune the cavity into resonance and shut 
off the external trap. 



In addition to investigating the dynamics of the atom 
in the presence of the standing wave mode, we have sug- 
gested various applications for state preparations and 
measurements. It has been explained and shown numer- 
ically how an incident atom may be reflected by the cav- 
ity while simultaneously leaving one photon into it, which 
might be used for creation of Fock states. The model can 
also be used for 'Stern-Gerlach' type of measurements, 
separating the upper and lower atomic states spatially. 
Likewise, it could also be used for separating the vacuum 
|0) from the one photon state |1) and to create various 
entangled states, like EPR-states. 

The paper is purely theoretical, and we understand 
that experimental measurement of the effective masses 
will pose several challenges, depending on the approach 
chosen. First of all, the nodes of a free-space standing 
wave have a tendency to drift, which is avoided in a nat- 
ural way by using an optical or a microwave cavity. Both 
regimes, however, bring new difficulties: A periodic field 
mode is needed in order to observe changes predicted 
in the propagation characteristics; present high-Q mi- 
crowave cavities have, unfortunately, dimensions of the 
order of wavelength. Hence the optical wave number is 
not well defined and the Floquet approach adapted in this 
paper is no longer applicable |4jj. In the optical regime, 
time scales of both the atomic spontaneous emission and 
the cavity decay are usually shorter than the interaction 
time, and both processes will strongly affect the atom- 
cavity system in a way that is not included in the present 
consideration. One way to go around these difficulties is 
to use extremely long lived transitions and configurations 
in which upper atomic levels are only slightly populated, 
cf. Figs. and EJ the cavity decay may be circumvented 
by using a driven cavity, hence rendering the light field 
effectively classical. Similar experiments are performed 
using classical laser fields made stable by a set of mirrors, 
see e.g. The same kind of setup could in principle 

be used for determination of the effective mass and one 
may then neglect decay of the field. The dynamics of 
the atom will then be the same when it interacts with a 
classical field as with a Fock state. 

The Floquet approach presented for the description of 
atomic motion in terms of an effective mass parameters is, 
in principle, independent of weather the two-level atom 
interacts with a classical light field or with a photonic 
Fock state of a quantised field. The present results may, 
however, not be realistically achievable in today's cav- 
ity QED systems. Our treatment does set the stage for 
the concept of an electromagnetic effective mass and we 
may only hope that novel experimental conditions will, 
one day, make possible the observation of cavity-induced 
mass modification. 
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